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ABSTRACT 


A sinusoidally time-varying pattern of the values of the Newton’s eonstant of gravi¬ 
tation G measured in Earth-based laboratories over the latest deeades has been reeently 
reported in the literature. We put to the test the hypothesis that the aforementioned har¬ 
monic variation may pertain G itself in a direct and independent way. We numerically 
integrated the ad-hoc modified equations of motion of the major bodies of the Solar 
System by finding that the orbits of the planets would be altered by an unacceptably 
larger amount in view of the present-day high accuracy astrometric measurements. In 
the case of Saturn, its geocentric right ascension a, declination d and range p would 
be affected up to 10^ - 10^ milliarcseconds and 10^ km, respectively; the present- 
day residuals of such observables are as little as about 4 milliarcseconds and 10“^ 
km, respectively. We analytically calculated the long-term orbital effects induced by 
the putative harmonic variation of G at hand finding non-zero rates of change for the 
semimajor axis a, the eccentricity e and the argument of pericenter to of a test particle. 
For the LAGEOS satellite, an orbital increase as large as 3.9 m yr“' is predicted, in 
contrast with the observed decay of -0.203 + 0.035 m yr“k An anomalous perihelion 
precession as large as 14 arcseconds per century is implied for Saturn, while latest 
observations constrain it down to the lO”"* arcseconds per century level. The rejection 
level provided by Mercury is of the same order of magnitude. 


Subject headings: gravitation-celestial mechanics-ephemerides 
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1. Introduction 


The Newton’s gravitation al eonstant G (lGillieslll997l: iMohr. Taylor & Newellll^l2r) . 
measured for the first timeQ by Cavendish ( 1798 ) at the end of the eighteenth eentur\@, is o ne of 
the fundamental p^ameters of Nature setting the magnitude of the gravitational interaetion dUzan 
2003 . 2 OO 9 L 201 ll: Chiba 2011 1. In either the Newtonian and the Einsteinian theories it is assumed 
that it does depend neither on spatial nor temporal eoordinates, being a truly universal eonstant. 

_ In the twentieth eentury, prominent seientists ( Milnel 19351 1937 : Dirael 1937 : Jordan 1937 . 

1939h . mainly on the basis of eosmologieal arguments, argued that, aetually, G may experienee 
slow time variations over the eons. Current researeh took over sueh a fas cinating idea, so that 
nowadays th e re are several t heore ti cal scenarios encompassing it; s e e, e.g.. Brans & Dickel (1961 ); 
Brans ( 1962 ): Wu & Wang ( 1986 ): Ivashchuk & MeTNikov ( 1988 ): Melnikov! ( 2002 . 2009). 


Recently, lAnderson et al.l (l2015ah showed tha t a set of measurements of G obtained 


over the years with different techniques (see, e.g.. 

(Gundlach & Merkowitz 2000: Ouinn et al. 

2001: Schlammineer et al. 20061: Fixler et a 

2007 

Lamnoresi et a 

2008:IParks & Fallei 2010 

Luo et al. 20091: Tu et al.l 2010: Ouinn et al. 

2013; 

Rosi et al. 2014: 

Schlammineer 2014)) in 

terrestrial laboratories over the latest decades (Speake & Ouinn 2014|) can be satisfactorily 


modeled with the following harmonic time-dependent signature 


with (lAnderson et al.ll2015al) 


) — Go -i- AG(t) — Go -i- Aq sin {uiQt (^) , 

(1) 

Go = 6.673899 x 10““ kg'^ m^ s"^ 

(2) 

Ag = 1.619 X 10“^^ kg-' m^ s-^ 

(3) 

^ = 80.9 deg. 

(4) 

2n 

(5) 

“ G 


^For a critical discussion of ea rlier attempts (iBougued 1 17491: iMaskelynd Il775h . see, e.g.. 


Poyntingl 1 8941) : ISmallwoodI (120 1 01) . 


^In view of the unit conventions in use at the time of Cavendish’s work, the gravitational con¬ 
stant did not appear explicitly in it. Indeed, one of the first explicit formulati ons of the New¬ 
tonian gravity in terms of a gravitational constant G appeared not before 1873 dCornu & Bailie 


1873 ). Thus, the English scientist actually measured the mean density of the Earth dClotfeltei 


1987 ) . His results were later refor mulated in modem terms as a determination of G dPoynting 


1894 : Mackenzie 1900l: Boys! 1893) . 
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Pg = 5.899 yr. 


( 6 ) 


Anderson et al. ( 2015al) did not suggest that the reported pattern may be due to some modifieations 
of the currently accepted laws of gravity. They remarked that a corr elation with recently repo rted 
harmonic variations in measurements of the Length of Day (LOD) (iHolme & de Vironll2013h is 
present. Later, a compilation of all published measurements of G performed since 1980, inclu ding 
also so me additions an d corrections with respect to the value s analy zed in lAnderson et al. 


atso so me aaartrons an a correctrons wrtn respect to me vatue s anaty zea rn lAnaerson et at. 
i2015ah . was offered by lSchlamminger. Gundlach & NewmanI 1201511 . Their least-square fit to 
such an expan ded dataset confirmed t he existence of the sinusoidal component with a 5.9 yr 
periodicity by lAnderson et al.l (l2015ah . although the correlation with the LOD turned out to 
be weakened; in the minimization procedure , a second pe riodicity of about 1 yr was found by 
Schla mminger. Gundlach & Newman ( 2015b . Cautiously, Schlamminger. Gundlach & Newman 


(1201 5b did not explicitly mention any possible variation of the gravitational laws, warning 
that the harmonic pattern they found r nay be due to some-unspecified-systematic errors, with 
un derestimated systematic uncertainty. Pitkin j[2015b . by using the compiled measurements of G 
by I Schlamminger, Gundlach & NewmanI (l2015b and a Bayesian model comparison, claimed that a 
constant G measurement model with an additional Gaussian noise term would fit the experimen tal 


data be tter than a model containing periodic terms. In the re ply they offer ed, lAnderson et al. 


( 2015b ) were not able to independently confirm the claim by Pitkinl ( 2015 ). and stood by thei r 
conclu sions of potential periodic terms in the reported G measurements. In particular. IPitkin 
( 2015b found that a mo del with two hamionic c omponents would be slightly favored with respect 
to their early proposal ( Anderson et al. 20L5ab . al lowing also to cope with the issue of the 


weakened correlation with the LOD pointed out in I Schlamminger. Gundlach & NewmanI (l2015l) . 


At the time of writing, none of the putative systematic errors allegedly affectin g the 


Earth-based measurements of G, reasonably conjectured by lAnderson et al.l (l2015al) and 
Schlamminger, Gundlach & Newman ( 2015 ). have been yet discloseci^. Things standing thus 
unknown, it is not unreasonable to follow a complementary approach, which could help in 
effectively selecting the directions of further experimental analyses, and look at the potentially 
intriguing-although admittedly unlikely-possibility that some modifications of the curre ntly 
accepted laws of gravity may be at work in the present case. According to KleinI ( 2015 ). the 
observed discrepancies between the G values determined by different experiments may be 
connected with a differential interpretation of MOND theory applied to the galaxy rotation curves. 


In this spirit, we want to quantitatively put to the test the daring conjecture that (some of) 
the observed sinusoidal pattem(s) in the G data records may reflect an unexpected physical 
phenomenon in an independent way by looking at the consequences that such an effect, if real, 
would have on systems other than those used to collect the measured values of G on the Earth. To 
this aim, for the sake of simplicity, we will consider the simpler model of eq. ([IJ-eq. (jb]) and the 


^A search for them is outside the scope of the present work. 
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changes which would occur in the motions of the major bodies of the Solar System to eheek if 
they are compatible with the eurrent stringent limits posed on their standard dynamics by accurate 
astrometric measurements. 


Here, we will use reeently released Cassini data analyses spanning ten years (20 04-2014) of 


the orbit of Saturn in te rms of its geocentrie range p, right aseension a, deelination 6 (IHees et al. 


20l4IJones et al.ll2015h to independently test the G(t) seenario of eq. Q-eq. db]). In partieular, we 


will suitably eompare numerically simulated signatures Aa(t), A6(t), Ap{t) induced by eg. 
dbl) on the Kronian Cele stial eoordinates with the eurrently existing residuals for them (IHees et al 


1^14 : Jones et aP 2015b . Full details of the methodology adopted are given in Seetion|2l Seetion 
|3] summarizes our findings. 


2. Methods and results 

A striking feature of the alleged time-variation of G investigated here is its relatively short 
oharaeteristie time seale which, according to eq. d^, amounts to just about 6 yr. This is in 
neat eontrast with virtually all the theoretieal models predicting a G{t) varying over typieally 
eosmologieal timeseales. This distinetive feature has also direet phenomenologieal consequenees. 



may not be straightforwardly extended to the present ease since they were inferred from 
least-square reduetions of planetary and lunar positional data by modeling AG(t) as a seeular 
trend. Sueh a ehoiee, reasonable in view of the extremely slow ehanges assumed in the literature 
for G with respeet to the typieal orbital frequeneies of the major bodies of the Solar System, does 
not apply to eq. dD- Thus, a dedicated analysis should be performed in the present ease: it will 
be the subjeet of the present Seetion. In particular, in Section [2T] a numerical approach will be 
followed, while an analytieal ealeulation will be offered in Seetion |2^ 


2.1. A numerical approach 

As a first step, we simultaneously integrate the baryeentrie equations of motion of all of the 
eurrently known major bodies of the Solar system in reetangular Cartesian eoordinates over a 
eentennial time-span (1914-2014) with the standard paekage MATHEMATICAL The dynamieal 


"^The MATHEMATICA method adopted is ExplieitRungeKutta, with a working preeision used 
in internal eomputations of 53 logjQ 2 « 16, and 8 digits of preeision and of absolute aeeuraey. The 
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accelerations modeled inelude the General Theory of Relativity to the first Post-Newtonian 
level, and all the major known Newtonian effeets like the Sun’s oblateness, pointlike mutual 
perturbations by the eight planets and the t hree largest asteroids, two mas sive rings aeeouriting f or 
the minor asteroids ( Kuehynka et al. 2010l) and the Kuiper Belt’s objeets ( Pitjeva & Pitjev 2013 1. 

The initial eonditions are taken from Tables [H to [2j they come from a n adjustment of the su ite 
of measurement and dynamical models of the EPM2013 ephemerides (IPitjeva & Pit) evl 120141 to 
an extended data record of more than 800 000 observations ranging covering last century, and are 
referred to the epoch to = JD 2446000.5 (27 October 1984, h: 00.00.00). 

Thus, keeping the other parameters of the numerical integration unchanged, we repeat 
the same step by including also the putative variation of G according to eq. Q-eq. db]). Both 
numerical integrations, with and without AG(t), share the same initial conditions for the known 
bodies of the Solar system retrieved from Tables [U to [2l From the resulting time series of the 
Earth and Saturn, we numerically compute the time series of p, a, 6 of Saturn, with and without 
AG(ty, then, for all of the three Kronian Celestial coordinates, we compute differential time 
series Aa(t), A6(t), Apit) which show up the impact o f AG{t) over the 2004-2014 interva l of 
time covered by the most recent Cassini data analyses (IHees et al.ll2014l:ljones et al.ll2015h : also 
Gaussian white noise is added to properly simulate the impact of the measurement errors. The 
results are displayed in Fig. [T] Finally , we compye our simulated resi duals to the corresp onding 
existing Kronian residuals in Fig. 5 of Hees et aP (2014) and Fig. 2 of Jones et al. ( 2015b . which 
were produced without explicitly modelling the perturbing action of AG(t). 

Some technical considerations about the general validity of the approach followed are 
in order. Possible objections of lacking of meaningfulness concerning such kind of direct 
comparisons among theoretically calculated signatures of a certain dynamical effect and actual 
data processed without modelling the effect itself have recently proved to be ineffective, at least in 
some specific cases. Indeed, apart fro m the fact that such an approach had been proven successfu l 
since the time of the Pioneer anomaly ( lorio & Giudice 2006 : Standish 2008 . 2010l: Fienga et al. 


2010h . the latest constraints on a certain form of the MOND theory, equivalent to the action 
of a remote trans-Plutonian body located in the direction of the Galactic C enter, which were 
obtained by explicitly modelling it in a dedicated planetary data reduction (IHees et al.ll2014h . 
turned out to be equivalent to those previously established by comparing theoreticall y cornputed 
effects to their observationally inferred counterparts determined without modelling it (Iloridl2010b . 
Finally, one might raise consisten cy issues about pur ari alysis since our initial condit ions come 


from the EPM2013 eph emerides ( Pitjeva & Pitjev 2014b . while the post-fit residuals ( Hees et al. 


20141: iJones et al.ll2015h to which our simulated signatures are contrasted were obtained with 


the DE430 ephemerides. Actually, it is not so. First, the EPM2013 initial planetary state 
vectors differ from the DE431 coordinates at to, retrieved from the HORIZONS Web interface at 


http://ssd.ipl.nasa 

.gov/?horizons, by just e = 2 (xrpm - Vur) / (vrpm + Vur) » 10 Then, from 

Fig. 7 of 

Pitieva 

(2013) it can be noticed that the differences between the Kronian Celestial 


computer used has a 64-bit operating system. 
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Table 1: Solar System Barycentric (SSB) initial positions jcq, yo, Zo of the Sun, the eight p lanets 
and the dwarf planet Pluto estimated with the EPM2013 ephemerides ( Pitjeva & Pitjevll2014h . They 
are referred to the International Celestial Referenee Frame (ICRF2) at the epoeh JD 2446000.5 (27 
October 1984, h: 00.00.00). (E.V. Pitjeva, private communication). 



xo (au) 

yo (au) 

Zo (au) 

Sun 

0.0005203216237770 

0.0084630932909853 

0.0035304787263118 

Mercury 

-0.1804160787007675 

-0.3766844130175976 

-0.1834161566766066 

Venus 

0.3542041788200310 

-0.5632690218672094 

-0.2760410432968001 

Earth 

0.8246466037849634 

0.5178999919107230 

0.2244221978752543 

Mars 

1.1810208746662094 

-0.6327726058206141 

-0.3225338798256297 

Jupiter 

1.6184148186514691 

-4.4994366329965958 

-1.9681878191442090 

Saturn 

-6.5110893457860746 

-6.9749292526618341 

-2.6006192842623910 

Uranus 

-5.3939091822416660 

-16.7458210494500221 

-7.2578236561504630 

Neptune 

0.5400549266512163 

-27.9860579257785673 

-11.4683596423229197 

Pluto 

-24.1120716069166541 

-17.4030370013573190 

1.8337998706332588 


Table 2: SSB initial velocities xq, yo, Zo of the major bodies of the Solar system. The other details 
are as in Table [B 


Sun 

Mercury 

Venus 

Earth 

Mars 

Jupiter 

Saturn 

Uranus 

Neptune 

Pluto 


io (au d 

-0.0000082057731502 

0.0202454878200986 

0.0175313054842932 

-0.0099025409715828 

0.0078002404125410 

0.0070681002624717 

0.0038982964656769 

0.0037423336376372 

0.0031184897780294 

0.0019439712118417 


yo (au d‘^) 

-0.0000014976919650 

-0.0077363989584885 

0.0093063903885975 

0.0130300000107108 

0.0120303517832832 

0.0025516512292944 

-0.0033491220519229 

-0.0011699461011685 

0.0000964171123061 

-0.0025828344901677 


Zo (au d“') 

-0.0000004352388239 

-0.0062333734254165 

0.0030763136627655 

0.0056499201171062 

0.0053067080995448 

0.0009215336853663 

-0.0015508256594232 

-0.0005654073338528 

-0.0000381527005270 

-0.0013916896197374 
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Saturn: right ascension 



Saturn: declination 



Saturn: range 



Fig. 1.— Numerically produced time series of p, a, 5 of Saturn perturbed by AG(t) as in eq. O- 
eq. dH). For each of the three Kronian observables considered, they were calculated as differences 
between two numerical integrations of the SSB barycentric equations of motion of the Sun, its eight 
planets and the dwarf planet Pluto from 1914 to 2014 with and without AG(0- Both integrations 
shared the same initial conditions in rectangular Cartesian coordinates, retrieved from Tables [H to 
|2l and the same standard dynamical models, apart from AG(t) itself. Thus, such curves represent 
the expected AG(0-induced signatures Ap, Aa, Ad. The pattern s and the size of the present signals 
can be c ompared with the r ange residuals bv iHees et al.l(l2014h (Apexp <0.1 km) and those for RA, 
DEC by Jones et aP ( 2015 ) (Aorexp, Adexp < 4 milliarcseconds). 
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coordinates calculated with the EPM2011 and DE424 ep hemerides oyer a t i me interval as little 
as 10 years (2004-2014) are smaller than the residuals in Hees et aP (2014); Jones et al. ( 2015 1 
obtained with the DE430 ephemerides. 


2.2. An analytical calculation 


The sinusoidal part of G induces a time-dependent component of the gravitational 
acceleration which, to the Newtonian level, can be written as 


( 8 ) 


In view of eq. dS]), eq. ([8]) can be thought as a small radial correction to the usual inverse- 
square law. As such, it can be treated vyith the s tandard methods of the perturbation theory 
( Bertotti. Earinella & Vokrouhlickil2003 : 3^12008 ) by assuming a Keplerian ellipse as unperturbed, 
reference trajectory. 

Eor a test particle orbiting a primary of mass M along a Keplerian ellipse with semimajor axis 
ao, orbital period Pq = Inri^^ = 2n M~'^ and eccentricity eo, the time t is connected with 

the tr ue anomaly /, which gives the instantaneous position along the orbit, through (ICapderou 
2005[) 


no 


{t = 2 arctan 


1-^0^ f 

-tan - 

1 + Co \2 


Co 




Cq sin/ 


(9) 


1 -I- Co cos / 

In eq. (|9l), hq is the unperturbed Keplerian mean motion, while tp is the time of passage at 
pericenter. The true anomaly / usually starts to be reckoned just at the crossing of the pericenter, 
so that / = 0. Eor the sake of simplicity, in the following we will assume tp = 0. 

By inserting eq. ([8]), with eq. dD and t given explicitly by eg. (|9l), in the right-hand sides of 
the Gauss equations (IBertotti. Earinella & Vokrouhlickll2003l:IXull2008l) for the variation of the 
Keplerian orbital element^ 

da 2 r I Pn\ 

coA«sin/ + Ar ^ , (10) 


dt 


Hq 




ro 


de 

dt 


v> 


tioao 


sin / -I- At 


Co \ ^0 


( 11 ) 


^In eq. (fT0l) -eq. (fT3l) . Ar, At, are the components of the disturbing acceleration A onto the 
radial, transverse and out-of-plane directions, Q. is the longitude of the ascending node, co is the 
argument of pericenter, and I is the orbital inclination to the reference {x,y} plane. 
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dt 


riQao 




el sin/o 


Aa, I — I sin (ojo + f ), 
ao 


( 12 ) 




d(o 

dt fiQGQeQ 


,2 , 


-Ar cos/ + |1 + ^|sin/ 


, dQ 
cos/o-—, 
dt 


evaluated onto the unperturbed Keplerian ellipse 


ro = 


Po 


Po 


= ao(l -el'j. 


1 + eo cos / 

and averaging them over one orbital period Pq by means of dCapderou 2005 ) 

v3/2 


dt 


(l - ‘■o) 


df no(l +eocos/)^’ 


(13) 


(14) 


(15) 


non-vanishing long-term rates of ehange of the semimajor axis, the eeeentricity and the perieenter 
are obtained. They are 


I da 

i dt 


eoAcMPoPc Jt^Pq 


jf-al 


[n - pi) 


cos 


-I- (/> sm 




(16) 


1 dt I 


ArMPnP. 


O^G 


An^al [Pl - Pi) 


sin0 


'InPo 

sm I-h (p 


+ O (^o), 


(17) 


\ dt I 


ArMPlP. 


O^G 


jp-al 


(pl - O 


sm 


sin (+ <p\ + 0 (eo) • 


(18) 


While the leading component of the rate of change of the semimajor axis is of order <9 (eo), the 
next-to-leading order terms of the rates of the eccentricity and the perieenter are of order (9 (eo)- 

To the approximation level in e indicated, which is quite adequate in the Solar System, eq. 
(fT^-eq. (fTSl) are exact in the sense that no a priori assumptions on the relative magnitudes of 
Pq and Pq were taken. Thus, they can be applied to a variety of systems ranging, e.g., from the 
fast Earth’s artificial satellites to the slowest planet of the Sun. As such, eq. (fT^-eq. (fTSl) hold 
also for any putative modified model of gravity yield in g possibly a harmonic var i ation of G ove r 


arbitrary timescales (iMorikawall 19901: IBarrow 


Barrow & CarnI 19961:iBarrow & Parsonslll997 


1993b a: Barrow & Mimoso 1994 : Barrowl 1995 


Stadnik & Flambaumll2015 ). 
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Let us start with the LAGEOS satellite orbiting the Earth along a nearly circular orbit 
with an = 12274 km, eo = 0.0039, Eq = 3.7 hr = 4.2 x lO'"* yr. It is known since long time 
( Rubincainl 1982 ) that its semimajor axis experiences a s ecular decrease due to a variety o f 
physical mechanisms. Its latest measurement amounts to ( Sosnica et al. 2014 : Sosnical2014h 


^LAGEOS 


= -0.203 + 0.035 m yr 


.-1 


(19) 


Instead, eq. (fT^ predicts a secular increase as large as 

=3.9myr-'. 


LAGEOS 


( 20 ) 


which is in disagreement with eq. (fT9l) . 

Moving to Saturn, which orbits the Sun in about 29 yr at 9.5 au, it turns out that eq. (fT^ and 
eq. (fTSl) predict 


. pred 
7 


-668 myr“\ 


( 21 ) 


= 14 " cty-\ 


( 22 ) 


Actually, there is no trace of such an enormous secular decrease of the Kronian semimajor axis in 
the observational data records. As far as its per ihelion is concerned, the current admitted range for 
a putative anomalous precession is as little a^ ( Pitjeva & Pitjev 2013 : Eienga et al. 2011 ) 


Adr^™ = (-3.2 + 4.7) x lO”'^ " cty-\ 


(23) 


Adr; 


INPOP 


1 ) 


= (1.5 + 6.5) X 10'^ " cty"'. 


(24) 


Similar results hold also for other planets like, e.g.. Mercury. His predicted perihelion rate is as 
large as 

=-108"cty-', (25) 


while any anomalous precession is constrained within (IPitjeva & Pitjevll2013l: lEienga et al.ll2011() 

Adr™ = (-2 + 3) X 10"^ " cty-\ (26) 


Adri 


• INPOP 


(4 + 6)x 10'^ " cty-'. 


(27) 


^Here, tzt = G + m is the longitude of pericenter. 
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3. Summary and conclusions 

In this paper, we have independently tested the hypothesis that the harmonic temporal 
variation in the time series of the laboratory measurements of the Newtonian constant of 
gravitation made over the years, which has been recently reported in the literature, may be due 
to some unknown physical mechanism affecting G itself. Cautiously, its discoverers did not 
mention such a possibility, claiming instead that it may be likely due to some systematic errors. 
Nonetheless, at present, none of them has been either identified or even explicitly suggested. 

To this aim, we looked at the effects that such a putative time-dependent behaviour of 
the fundamental parameter characterizing the strength of the gravitational interaction, if real, 
would have on the orbital motions of the major bodies of the Solar System. We numerically 
integrated their equations of motion with and without the proposed modification, and calculated 
the differences of the resulting time series for the observables used in real astrometric data 
reductions (right ascension a, declination d, range p) to produce simulated residual signals 
Aa(t), AS(t), Ap{t). We remark that, given the specific functional dependence of the effect 
considered and its relatively short characteristic time scale compared to the typical Solar System’s 
orbital periods, it would be incorrect to straightforwardly extend the existing bounds on G/G to 
the present case since they were obtained by modeling a secular variation of G. It turned out that 
the resulting anomalous signatures in right ascension, declination and range are far too large to 
have escaped detection in the residuals produced so far with the existing standard ephemerides, 
even if the putative variation of G was not explicitly modeled in all of them. Suffice it to say that, 
in the case of Saturn, cr, 6 and p would be affected up to 10"^ - 10^ milliarcseconds and 10^ km, 
while the current residuals are as little as about 4 milliarcseconds and 10“* km, respectively. 

Complementarily, we performed also an analytical calculation of the long-term perturbations 
which a harmonic variation of G would induce on the orbital elements of a test particle orbiting 
a central mass. In the limit of small eccentricities, we obtained non-vanishing secular variations 
of the semimajor axis a, the eccentricity e and the pericenter oj. They have a general validity 
since no a priori assumptions concerning the relative magnitudes of the orbital and G frequencies 
were made; as such, they can be applied, in principle, to any physical mechanism predicting 
a sinusoidally time-dependent variation of G. A comparison with the latest observational 
determinations for the LAGEOS satellite and Saturn confirmed the outcome of the numerical 
analysis. Indeed, the recently measured orbital decay of the geodetic satellite, of the order of 
-0.2 + 0.03 m yr“\ rules out the predicted increase of « -1-4 m yr“^ Furthermore, the expected 
anomalous apsidal rate of 14 arcseconds per century for Saturn falls neatly outside of the allowed 
range for any possible anomalous Kronian perihelion precession, which is set by the observations 
to the « lO"'* " cty“' level. Also Mercury yields an analogous outcome since the magnitude of 
its predicted anomalous perihelion precession amounts to 108 " cty“\ while the observations 
constrain any possible deviations from standard physics down to « 10~^ - lO-"* " cty“^ 

In conclusion, our analysis quantitatively rules out the possibility that some long-range 
modification of the currently accepted laws of the gravitational interaction may be at work in the 
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present ease aeeounting for the observed harmonie pattern of the laboratory-measured values of G. 
As such, it may contribute to further direct future investigations towards the discovery of the even 
more likely systematic uncertainties allegedly plaguing the set of measurements of G analyzed so 
far. 
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